Origin of Long Lived Coherences in Light-Harvesting Complexes 
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A vibronic exciton model is developed to investigate the origin of long lived coherences in light- 
harvesting complexes. Using experimentally determined parameters and uncorrelated site energy 
fluctuations, the model predicts oscillations in the nonlinear spectra of the Fenna-Matthews-Olson 
(FMO) complex with a dephasing time of 1.3 ps at 77 K. These oscillations correspond to the coher- 
ent superposition of vibronic exciton states with dominant contributions from vibrational excitations 
on the same pigment. Purely electronic coherences are found to decay on a 200 fs timescale. 



The role of quantum mechanics in biological processes, 
as well as photosynthetic light-harvesting, has been of 
interest for a long time The topic received re- 

newed attention after the observation of long lived oscil- 
lations in the two-dimensional (2D) spectra of the Fenna- 
Matthews-Olson (FMO) protein pigment complex 0]. 
These oscillations were interpreted as a signature of elec- 
tronic coherences between the delocalized energy eigen- 
states of the complex, and it was argued that their slow 
dephasing could enhance the efficiency of energy transfer 
between the chlorosome antenna and the reaction center 
0, [j| . Subsequent studies revealed that the oscillations in 
FMO have a dephasing time of 1.2 ps at 77 K [6(, and that 
such oscillations are a common feature in light-harvesting 
complexes 0, Q . Based on the known structure of FMO 
(Fig. [TJ, simulations employing formally exact equations 
of motions for the reduced density operator found oscilla- 
tions with dephasing times of 200 — 300 fs - clearly shorter 
than the experimental observation [El [9l-[Tl| . Significantly 
longer dephasing times were found if the transition fre- 
quency fluctuations (static or dynamic) of the different 
pigments are assumed to be correlated [12, EH- However, 
theoretical studies using molecular dynamics simulations 
of the interaction between electronic and nuclear degrees 
of freedom (DOF) have not been able to confirm the ex- 
istence of such correlations [nt or long lived electronic 
coherences [IE[- To what extent the long lived oscilla- 
tions in the experiments reflect electronic coherences, or 
if they influence the transport of energy across the com- 
plex 0,113, 

thus remains an open question. 
Analysis of excitation dynamics in molecular aggre- 
gates typically employs a reduced description, where the 
electronic DOF and their mutual couplings are treated 
explicitly, and the nuclear modes of the pigments and 
protein are treated as a heat bath [13, EH. The pres- 
ence of underdamped vibrational modes in the bath pro- 
duces oscillatory signatures in 2D spectra which are simi- 
lar to the modulations predicted for electronic coherences 
[l9j . The dephasing time of such nuclear coherence are 
of the order of several picoseconds, and they can often 
be treated as completely undamped on the time scale 
of a typical 2D experiment. Low temperature fluores- 
cence line narrowing (FLN) experiments on FMO have 



revealed a large number of vibrational modes in the range 
of — 350 cm _1 [20j. However, the oscillations seen in the 
2D spectra of FMO cannot be directly related to simple 
nuclear wavepackets, because the frequency of the oscilla- 
tions does not match any of the vibrational frequencies, 
and the Huang-Rhys factors of the modes are too low. 
On the other hand, recent simulations have revealed un- 
expected effects on the electronic structure and dynamics 
if vibrational modes are explicitly included in the system 
[2ll | . Motivated by these results, we develop a vibronic 
exciton Hamiltonian in so called one particle approxima- 
tion [22I, [23| for FMO, in which one vibrational mode on 
each monomer is treated explicitly. We will show below 
that this model predicts oscillations in the 2D spectra of 
FMO with 1.3 ps dephasing times at 77 K, and that these 
long lived coherences can be traced to superpositions of 
vibronic exciton states located on the same pigment. 

The total Hamiltonian of a molecular aggregate in con- 
tact with the environment is partitioned in a standard 
way into system, bath and system-bath interaction terms, 
H = Hs + Hb + Hs b ■ The system Hamiltonian describes 
the Q y transition on each bacteriochlorophyll (BChl) in 
FMO (see Fig. [1} with the vibrational progression of a 
single vibrational mode. Including the resonance cou- 
pling between the transitions, the system Hamiltonian 
reads 
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where E n is the transition frequency of pigment n (site 
energy), wo is the vibrational frequency, and v is the 
quantum number of the vibrational mode. The cou- 
pling energy J n ,i/;m,u' between the individual transitions 
can be expressed via the electronic resonance coupling 
Jnm an d the Franck-Condon amplitudes of the vi- 
brational mode [1^| J n ,v,m.v' = ( i/ |0)(i / '|0). The eigen- 
values and wave-functions of Hs are given by Huj a and 
\ a ) = J2 n v c n.J n ' v )i respectively. The bath Hamilto- 
nian, Hb , is described as a collection of independent har- 
monic oscillators, for which the system-bath interaction 
is given by H S b = J2 n ,u ud„(uj)q n \n, v)(n, v\. Here q is 
a generalized coordinate of the environment, and d(u>) 
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is the displacement of the excited state relative to the 
ground state. We assume that the system-bath Hamilto- 
nian does not depend on the state of the vibrational mode 
(v), implying that a vibrational coherence on an isolated 
monomer is undamped. Assuming equal but uncorre- 
cted system-bath interaction for the different pigments, 
the energy gap correlation function in the local basis is 
given by C nm (t) = Lo 2 d n d m (q n (t)q m (0)) = C {t)8 nm . 
When the interaction of the system with the environ- 
mental modes (i.e. all except those treated explicitly in 
the system Hamiltonian) is weak, it is advantageous to 
perform calculations in the eigenstate basis of the system 
Hamiltonian. The correlation function of the energy gap 
in the eigenstate representation can be expressed via the 
expansion coefficients, „, and the correlation function 
of each excitation in the local basis Co(t), 

c^W = c (t){E„,,«,J 2 (O 2 

+ E [«,,) 2 (<^) 2 + « x ) 2 «J 2 ]}=CoW 7 ^ 

The dephasing dynamics of a coherent superposition be- 
tween the eigenstates a and j3 is determined by the line- 
shape function g a p(t) — f dr JZ dr / 7 Q , ( gCo(T / ). The cor- 
relation function Co(t) is connected to the the spectral 
density C"(uj) via a Fourier transform (see e.g. Ref. 
[2^|). In addition to dephasing, the system- bath inter- 
action leads to relaxation between the eigenstates of the 
system. We use Redfield theory [l?} (Markov and secular 
approximation), where the relaxation rate is given by 

x{(l + n(u a 0))C''(u a p)+n(-u) a p)C"(-u) a p)}, 

where C"(uj a p) > and n{uj a p) is the Bose-Einstein dis- 
tribution function. The total relaxation rate from a level 
a, T a = ^E 7 ^ a ^a->7' determines the lifetime broad- 
ening of this exciton state. Assuming k^T < hug, the 
linear absorption spectrum can be calculated as [18| 
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where a runs over all exciton levels and the tran- 
sition dipole moments fi a o are given by fi a o = 
Hn.v c n,v^n{v\0)- Here (. . .)a,o denotes the average over 
a random distribution of pigment energies and orienta- 
tions of complexes. 

To simulate the oscillations in a third order experiment 
(i.e. 2D spectra) we adopt the doorway window repre- 
sentation [25j. Of all Liouville pathways contributing to 
the signal, only those involving a coherence between two 
levels in the excited state will give rise to oscillations 
during the waiting time t 2 . Without the loss of gener- 
ality, we focus on the non-rephasing coherence pathways 
illustrated in Fig. [TJ which give rise to oscillations along 



the diagonal in the non-rephasing 2D spectrum The 
response function for this pathway is given by 

Ra0,O = (((/X Q 0) 2 (M/30) 2 )oG , Q (t 3 )G , 2](t2)G Q (tl))A, (5) 

where G a (t) = e -^ot-r a t-g aa (t) ; and G^(f 2 ) = 
e -iu a (,t2 e -g a c*{t2)-gi3i3{t2)+2g a fi(t2)~(.T a +r l j)t2 _ j n tmg wor ]j 

we use the site energies (E n ) and resonance couplings 
(Jnm) for FMO Chlorobium tepidum from Ref. [18| . and 
the analytical formula for the overdamped part of the 
spectral density, C"{uS), extracted from a FLN experi- 
ment (20| . The direction of the transition dipole moments 
were taken from the protein data bank file 3ENI (26| . 
FLN experiments have identified 30 vibrational modes in 
FMO, and the strongest feature in the_ spectrum arises 
from three modes around 185 cm -1 [2(3]. To retain a 
simple description, we treat this cluster of modes as one 
effective mode with a frequency of ojq = 185 cm -1 and 
a Huang-Rhys factor of 0.05. For calculations with the 
exciton model, the vibrational mode was included as a 
underdamped contribution in the spectral density. In 
all calculations presented in this paper, we sampled the 
pigment transition energies from a Gaussian distribution 
with a FWHM of 80 cm- 1 . 

The vibronic exciton and exciton model predict very 
similar linear optical properties as illustrated by the sim- 
ulated linear absorption spectra shown in Fig. [TJ Figure 
[2ta) shows the time evolution of coherences involving the 
lowest state of the exciton model. During t 2 the signals 
oscillate with frequencies corresponding to the splitting 
between the exciton levels and are completely damped af- 
ter 400 fs. The strong damping can be readily understood 
from Eqs. and For the exciton model we find that 
the cross-correlation term g a p is small and the coherence 
pathways decay mainly with exp(— g aa {t2) — gppfo))- If 
a (static or dynamic) correlation of the transition fre- 
quency fluctuations in the site basis is assumed (ad hoc) , 
the cross-correlation functions in the exciton basis be- 
come larger and enable longer dephasing times. 

The coherence pathways involving the lowest levels of 
the vibronic exciton model are shown in Fig. [^Jb). The 
vibronic exciton coherences are remarkably long lived, 
and the signals show only minor damping on a 2 ps 
timescale. The long dephasing time of the coherences 
in the vibronic exciton model can be understood by in- 
spection of the expansion factors (|c" „| 2 )a given in Tab. 
UJand the transition frequency distributions shown in Fig. 
[TJ For instance, state 1 corresponds to 75 % to an excita- 
tion of the v = transition of pigment 3, while state 4 has 
a large contribution of vibrational excitation [y — 1) on 
the same pigment. As discussed above, the system-bath 
interaction is independent of the state of the vibrational 
mode, and these two vibronic exciton levels will there- 
fore experience highly correlated fluctuations, resulting 
in slow dephasing of the 1, 4 coherence. Despite the large 
contribution from the vibrational excitation to state 4, 
it has a transition dipole moment which is comparable 
to that of the other vibronic exciton levels. For non- 
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Figure 1: Linear absorption spectrum at 77 K for the exci- 
ton model (black) and the vibronic exciton model (red dash). 
The gray lines show the distribution of renormalized transi- 
tion frequencies weighted by the transition strength (scaled 
by 2/3) for state 1 and 4 in the vibronic exciton model. The 
filled areas illustrate the relative contribution from electronic 
(red) and vibronic (blue) excitations on pigment 3. The spec- 
tra of the vibronic exciton model have been shifted by the 
reorganization energy of the vibrational mode (—9.25 cm" 1 ) 
for comparison. The insets show the Feynman diagram il- 
lustrating the non-rephasing excited state coherence pathway 
in Eq. (12), and the arrangement of the 7 BChls in FMO 
(C. tepidum) [2|§|. This figure was generated using the VMD 
software (27l |. 



interacting pigments, only the zero-phonon state (v = 0) 
has a significant transition dipole moment. The strong 
transition dipole moment here is the result of intensity 
borrowing from the electronic transitions on the other 
pigments. As illustrated in Fig. [TJ state 4 of the vi- 
bronic exciton model exhibits both properties needed for 
the generation of long lived coherences: a large contribu- 
tion of vibrational excitation of pigment 3, and a strong 
transition dipole moment enabled by intensity borrowing. 
The combination of both effects leads to long lived oscil- 
lations on the red edge of the linear absorption spectrum 
as shown in Fig. [He). The oscillations show a bi-phasic 
behavior, where the initial 200 fs decay of the oscillation 
is due to the decay of coherences between vibronic ex- 
citon states localized on different pigments (like the 1,2 
coherence, see Fig. HJb) and Tab. HJ, while the long lived 
oscillations reflect coherences between vibronic exciton 
states localized on the same pigment. A fit to the oscil- 
lations give a dephasing time of 1.3 ps. By comparing 
to the stimulated emission signal calculated for the same 
spectral range, we estimate the modulation of the total 
signal, including stimulated emission and ground state 
bleach, to be 5 — 10 % for £2 > 0.3 ps. 

The Fourier transform of the signal in Fig. [2{c) is 
shown in Fig. HJd). The oscillation frequency of 205 
cm" 1 is higher than wo and also higher than the fre- 
quency observed in the experiment (160 cm" 1 ) [6]. The 
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Figure 2: Amplitude of the real part of the non-rephasing 
coherence pathways involving the lowest state, a) Re(Rip)n,A 
with ft = 2 (blue dash), /3 = 3 (red thin solid), and /3 = 4 
(black solid) for the exciton model, b) Re{R\f))n,& with f3 — 2 
(blue dash), /3 = 3 (red thin solid), and /3 = 4 (black solid) 
for the vibronic exciton model, c) Sum of all non-rephasing 
coherence pathways, ^2 a a Re(R a /3)sn,A , giving rise to signal 
in the range 12100 ± 30 cm -1 for the vibronic exciton model 
with vo = 185 cm -1 (blue) and t>o=117 cm -1 (red) at 77 K. 
The initial value of the signal is 0.48 and the first minimum at 
t2 — 0.04 ps has an amplitude of —0.23. The non-rephasing 
stimulated emission signal calculated for the same parameters 
and spectral range has an initial value of 0.17. d) Power 
spectrum of the Fourier transform of the signals in (c) starting 
from 0.2 ps. 



oscillation frequency depends on transition energies, elec- 
tronic couplings and vibrational frequencies according to 
Eq. JTJ). Figs, [^c) and (d) compares the oscillations on 
the red edge of the spectrum for two different vibrational 
frequencies found in the FLN experiment, ujq = 185 cm -1 
and wo = 117 cm" 1 . For ujq — 117 cm -1 , the oscillations 
have a frequency of 140 cm -1 and a shorter dephasing 
time as compared to ujq — 185 cm" 1 . The dephasing 
time of the coherences depend on the amount of vibra- 
tional character of the vibronic exciton states, and de- 
tailed analysis of the oscillations provide information on 
the energies and nature of the eigenstates not accessible 
from linear spectra (see Fig. [T]). However, the results in 
Figs. [2] (c) and (d) indicate that more than one vibra- 
tional mode needs to be included explicitly for a quanti- 
tative analysis of the oscillations in FMO. 

In this work we have shown that the vibronic exciton 
model predicts coherences in FMO with 1.3 ps dephasing 
times at 77 K. Our model does not invoke static nor dy- 
namic correlations in the site energies of the pigments, 
and uses experimentally determined spectral densities 
and vibrational frequencies. The long lived coherences 
are found to reflect coherent superpositions of vibronic 
exciton states with dominant contributions from vibra- 
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n = 3, v = 
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0.2 


0.03 
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0.67 


n = 4, y = 


0.21 


0.51 


0.01 


0.01 


n = 4, f = 1 


0.0 


0.0 


0.0 


0.01 


(/4}a(mIcm) 


0.87 


0.58 


1.3 


0.57 


(-B«)a - {-EiMcm" 1 ) 


0.0 


105 


175 


217 



Table I: Contributions of selected basis excitations (|c™ „| 2 ) 
to the four first vibronic exciton states averaged over ener- 
getic disorder. The numbering of the pig ments is defined in 
Fig. [T] and is the same as in Ref. [la ]. The two bottom 
rows show the averaged transition strength in units of the 
BChl monomer, and the average energy differences between 
the vibronic exciton levels, respectively. 



the system as well as on certain bath modes (Eq. (JTJ), 
both types of DOF become mixed. This enhances the ef- 
fective system-bath coupling in a way not accounted for 
in our exciton model. A similar mixing of system and 
bath DOF takes place implicitly when the reduced equa- 
tion of motion for the electronic system is propagated 
exactly [l^. It is clear that the mixing takes place for 
all bath modes. However, modes of the protein environ- 
ment are strongly damped and cannot contribute to a 
long dephasing time. As shown here, one or more of the 
underdamped vibrations found in the BChl monomers 
are needed to account for the enhanced dephasing times 
of the coherences. Our results imply that the oscillations 
in the 2D experiment on FMO reflect the dynamics of the 
nuclear DOF within a single pigment, which should have 
little impact on the transfer of energy from the chloro- 
some to the reaction centre. 



tional excitations on the same pigment. Because vibra- 
tional modes are an inherent property of all pigments, we 
expect vibronic excitons to be a general feature in the dy- 
namics of molecular aggregates. In the exciton language, 
the long lived coherences reported here correspond to a 
coherence between the system and the bath. Because the 
resonance coupling in the vibronic exciton model acts on 
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